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Abstract. We develop a systematic analytical approach on linear and nonlinear 
pulse propagations in an open A-type molecular system with Doppler broadening. 
O ' In linear case, by using residue theorem and a spectrum decomposition method, we 

CZ3 , prove that there exists a crossover from electromagnetically induced transparency 

(EIT) to Autler-Towncs splitting (ATS) for co-propagating configuration of probe and 
control fields. However, there is no EIT and hence no EIT-ATS crossover for counter- 
^ r-j propagating configuration. We give various explicit formulas, including probe-field 

Oh' spectrum decomposition, EIT condition, width of EIT transparency window, as well 

as a comparison with the result of cold molecules. Our analytical result agrees well 
with the experimental one reported recently by A. Lazoudis et al. [Phys. Rev. A 82, 
023812 (2010)]. In nonlinear case, by using the method of multiple-scales, we derive a 
nonlinear envelope equation for probe-field propagation. We show that stable ultraslow 
\Q ' solitons can be realized in the open molecular system. 
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1. INTRODUCTION 

In recent years, much attention has been paid to the study of quantum coherent 
phenomena in various multi-level systems, typical examples include Auter-Townes 
splitting (ATS) [I] and electromagnetically induced transparency (EIT) (2j. Such 
phenomena are not only important from viewpoint of basic research, but also very 
attractive for many practical applications, such as lasing without inversion, coherent 
population transfer, enhanced Kerr nonlinearity, slow light, quantum memory, atom 
and/or photon entanglement, precision spectroscopy, precision measurement, and so 
on [2113]. 

ATS occurs when absorption spectrum of a quantum transition can be decomposed 
into a sum of two net Lorentzian terms if one of two levels involved in the transition 
is coupled to a third level induced by a strong control field. EIT occurs when the 
absorption spectrum can be decomposed not only into two Lorentzians, but also with 
additional quantum destructive interference term(s). Usually, in systems with ATS or 
EIT, a transparency window is opened. However, the opening of the transparency 
window cannot be tell us whether the phenomenon belongs to ATS or EIT, each of 
which has different physical origin. ATS happens only for strong control field, but EIT 
happens even the control field is weak. Especially, Only for weak control field can 
essential characters of EIT be illustrated clearly [U El El tl] • 

EIT in various atomic systems has been studied intensively both theoretically and 
experimentally [21 [3]. However, systematic investigations of EIT in molecular systems 
are still lacking. Up to now, there are only several related experimental studies in 
molecular systems, including the EIT in 7 Li 2 [Sj, K 2 [10] and Na 2 vapors [TT1 [12], in 
acetylene molecules filled in hollow-core photonic crystal fibers [T3l E] and in photonic 
microcells [15], and in Cs 2 in a vapor cell [16], and so on. Major difficulties for observing 
EIT in molecules are small transition-dipole-moment matrix elements in comparison 
with those in atoms, and many decay pathways to other molecular states not involved 
in the main excitation scheme. 

In an interesting work reported recently by Lazoudis et al. [9], EIT in an open 
hot A-type molecular 7 Li 2 system has been studied experimentally. A numerical 
simulation under steady-state approximation is used by the authors for solving density 
matrix equations for molecules. Though the numerical simulation is helpful to explain 
experimental data, it is however hard to discern ATS from EIT objectively because the 
physical mechanism behind numerical results are not clear. In particular, since open 
molecular systems with Doppler broadening are very complicated and have very different 
features in different parameter regions, it is necessary to clarify in an analytical way the 
quantum interference characters inherent in such systems, which, to the best of our 
knowledge, has not been done in literature up to now. In addition, it is also necessary 
to go beyond steady-state approximation if probe pulse is used in experiment. 

In this work, we develop a systematic analytical approach on linear and nonlinear 
pulse propagations in open A-type molecular systems with Doppler broadening. In 



Analytical approach on linear and nonlinear pulse 



3 



linear case, by using residue theorem and spectrum decomposition method, we prove 
clearly that there exists a crossover EIT to ATS for co-propagating configuration of 
probe and control fields. However, there is no EIT and hence no EIT-ATS crossover 
for counter-propagating configuration. We provide various explicit formulas, including 
probe-spectrum decomposition, EIT condition, and width of EIT transparency window, 
as well as a comparison with the result of cold molecules. Our analytical result agrees 
well with the experimental one reported recently by A. Lazoudis et al. [9] . In nonlinear 
case, by using a standard method of multiple-scales, we derive a nonlinear envelope 
equation for probe-field propagation. We show that a stable ultraslow solitons can 
be realized in the open molecular system. Notice that nonlinear pulse propagation 
in coherent atomic systems via EIT has attracted tremendous attention in recent 
years [TT(, HHJ UHl EUJ ETJ |22} I23J [24] , nobody however has considered similar problem for 
molecules till now. 

The article is arranged as follows. In the next section we present our model and 
associated Maxwell-Bloch (MB) equations. In section |3l we consider the linear property 
of the system by using residue theorem and spectrum decomposition method. Quantum 
interference characters for hot molecules with both co- and counter-propagating 
configurations and also for cold molecules are analyzed in detail. In section HJ the 
method of multiple-scales is used to study the weak nonlinear propagation of the probe 
field. Lastly, section [5] contains a summary of the main results obtained in our work. 

2. Model 

The model adopted here is the same as that used in [9]. An open three-state A-type 
Li 2 molecular system (figured]) consists of an exited upper-level A 1 S+(w / = 5, J' = 13) 
(labeled |3)) and two ground states X 1 £+(?/' = 1, J" = 14) (labeled |1)) and X^+K = 
0, J" = 12) (labeled |2)). A control field with center frequency u c couples to the excited 
state |3) and the ground state |2). The other ground state |1) couples to the |3) by a 
probe field with center frequency u p . The exited level |3) decay spontaneously to the 
ground states |1) and |2) with decay rates and T 2 3, respectively. The parameter 7 
represents the transient relaxation rate of the molecule entering and leaving interaction 
region between light and the molecule. It reflects also the additional relaxation of each 
state due to the interaction with thermal reservoir [9]. The electric field vector of the 
system is E = ^2i= P c e i^i( z ^ i)e l ( ki ' r_£J(< )-|-c.c., where (k^) is the unit polarization vector 
(wave number) of the electric field component with the envelope Si (I = p,c). 

As indicated in the last section, decay processes in molecular systems are very 
complicated in comparison with those of atoms. There exist many decay pathways 
to other molecular states not involved in the main excitation scheme, and hence the 
theoretical model considered is necessarily an open one. In the excitation scheme 
adopted above, molecules occupying the excited level |3) may follow various relaxation 
pathways and decay to many lower vibration- rotation levels besides the levels |1) and 
|2). In our modeling all these levels are represented by the level |4). The decay rate 1^43 
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Figure 1. (Color online) A- type EIT scheme for open L12 molecular system. Excited 
state A^+O' = 5, J' = 13) (labeled |3)) couples to ground state X 1 ^^" = 0, J" = 
12) (labeled |2)) by the control field with center frequency uj c and also to another 
ground state X 1 E+(u" = 1, J" = 14) (labeled |1)) by the probe field with center 
frequency u p . A2 and A3 are detunings, Tji are population decay rates from \l) to \j), 
and 7 is transit rate. Molecules occupying the excited state |3) may decay to many 
other states besides the states |1) and |2). All these other states are represented by 
state |4). 



indicates the spontaneous emission rate of level |3) to level |4) (see figured]). 

For hot molecules, inhomogeneous Doppler broadening must be taken into account 
because the experiments are carried out in a heat-pipe oven [9] . The Hamiltonian of the 
system in interaction picture under electric-dipole and rotating-wave approximations is 

H = -h{Q c e ilkc - {r+vt) -^ t] \3}(2\ + fi p e i[kp - (r+vt) ^ 1 |3)(l| + c.c), (1) 

where v is molecular velocity, Q c (p) — ( e c( P ) • A*32(3i)) ^c(p)/(2^) is half Rabi frequency 
of the control (probe) field, with /x 7 the electric-dipole matrix element associated with 
the transition from state \j) to state |/). The optical Bloch equation in the interaction 
picture reads 

d 

+ ~ - ^13^33 + ^31 - ^pO"31 = °> 



^°"22 + ili.o-22 - 0-22) ~ ^23033 + Q* c o- 32 - Vl c al 2 = 0, 



d_ 

dt 
d 

^°"33 + «7(0-33 - ^33) + »T 3 «733 + ftpffgj + fi c CT* 2 

- n* p a 31 - tt* c o- 32 = 0, 
d 

i-Qi?4A + *7(°44 - C44) - ^43^33 = 0, 

d \ 

+ d 2 \ J <?21 + ^ c 031 - ^0-32 = 0, 
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d \ 

+ ^31 J 0"31 + ^p(o~ll - 0"33) + ^21 = 0, 

(i^ + d 32 ^j o- 32 + £l c (a 22 - 0-33) + Vt p a* 21 = 0, (2) 

for nondiagonal elements, where d 2 ± = — (k p — k c )-v+A2 — A1-H721, d 3 i = — k p -v+A3 — 
Ai+iTai, d 32 = -k c -v+A 3 -A 2 +z 732 with Ti , = (r j +r i )/2+ 7 + 7j c ° 1 (j,l = 1,2,3). Here 
Aj (j = 1, 2, 3) are detunings, and Fj denotes the total decay rate of population out of 
level \j), which is defined by Tj = ^ij- The quantity jjf is the dephasing rate due 
to processes such as elastic collisions. a^J is the thermal equilibrium value of o~jj when 
all electric-fields are absent. Equation (j5J) satisfies Ylij=\ a jj = 1 with X)j=i a j] = 1- At 
thermal equilibrium, population in the excited state |3) is much smaller than that of 
the ground states, i.e. (T33 ~ and hence <r^ + a 22 + C44 — 1- 

The evolution of the electric field is governed by the Maxwell equation. Due 
to the Doppler effect, the electric polarization intensity of the system is given by 
P = M a J^dvf(v){^ 13 a 31 exp[i(k p z-u p t)] + ^ 23 a 32 exp[i(k c z-uj c t)] + c.c.}, where Af a is 
molecular density and f(v) is the molecular velocity distribution function. For simplicity, 
we have assumed electric-field wavevectors are along ^-direction, i.e. k p c = (0,0, k pc ). 
Under the slowly-varying envelope approximation, the Maxwell equation reduces into 

( d 1 d \ f°° 
i \dz + cdi ) ^ p + Kn J dv f( v ) a 3i(z,v,t) = °> ( 3 ) 

with k 13 = Af a uj p \n 31 \ 2 /(2he c), here c is the light speed in vacuum. 

The MB equations (j2J) and ([3]) are our starting point for the study of linear and 
nonlinear pulse propagations in the open molecular system with Doppler broadening. 



3. Linear propagation 

3.1. Base state and general linear solution 



We first consider linear propagation of the probe field. For this aim, one must know 
the base state i.e. the 
Q p = 0. It is easy to obtain 



the base state i.e. the steady-state solution of the MB equations (J2J) and (j3J) for 



(0) _ [-fT^x, + (27 + r 43 )|ft c | 2 K 1 + r 13 |fi c | 2 (i - 

0~ 1 



'11 



X, 



(o) jjr^x, + \n c \ y% 

G n 



'22 



x 2 

(0) _ l\^c\-0- 2 i 



a 33 



CT (0) 
O44 



7|OjVg e q 
X 2 ' 
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and a { £ = = 0, where r 37 = 7 + T 3 , X 1 = {[(A 3 - A 2 ) - k c v} 2 + i 2 32 }/(2 l32 ) 
and X 2 = 7(7 + r 3 )Xi + (27 + r i3 + r 43 )|f2 c | 2 . Note that in above expressions 
d 2l = d 21 (v) = -(k p - k c )v + A 2 - A 1 + 2721, d 3i = d 31 (v) = -k p v + A 3 - A 1 + z'7 31) 
and d 32 = <i 3 2(f) = -k c v + A 3 - A 2 + «7 32 . 

When switching on the probe field, the base state (jl]) will be modified. In linear 
theory, Q p is taken as a very small quantity. At first order in Q p , the populations and 
the coherence between the states |2) and |3) are not changed, but with 



'8 



|fi c | 2 - (w + rf 21 )(u; + rf 31 ) 



— a 21 -F e , 
(i) = fa + d 2i)(o-n ~ o- { 3 f) + ^ c ^32 0) p ^ 

31 |(7 c |2_ (cj + d2l)(a; + rf3l) 

= a£ ) Fe» (5) 

where F is a constant, 9 = K(u)z — oot. The linear dispersion relation K(u) [25] is given 
by 

^ r * *, ^ fa + ^o^i? - 4?) + ^ 3 l 0) 



c J-^ \IL C \ 2 - (u + d 21 )(u + d 31 ) 

In thermal equilibrium, /(w) is the Maxwellian velocity distribution function, i.e. 
f(v) = l/(v / 7r vt) exp [— v 2 /v^}, with vt = \p2k^TjM the most probable speed at 
temperature T, and M the molecular mass. The integration in equation fl6]) with 
the Maxwellian distribution leads however to some complicated combination of error 
functions [26], which is very inconvenient for a simple and clear analytical approach. 
As did by Lee et al. [27], in the following we use the modified Lorentzian velocity 
distribution f{v) = vt l\\p^{v\ + v 2 )} to replace the Maxwellian distribution. 

We are interested in two different cases: co-propagating configuration (k p ~ k c ) 
and counter-propagating configuration (k p ~ — k c ), discussed below separately. 

3.2. Hot molecules with co-propagating configuration 

In this configuration, one has d 2 \ = A 2 — Ai + ij 2 i, d 3 i = —k p v + A 3 — A x + ?7 31 and 
^32 = — k p v + A 3 — A 2 + Z7 32 . The second term on the right-hand side of equation ([6]) 
can be calculated by using residue theorem [28]. There are two poles in the lower half 
complex plane 

k p v = A 3 — iX 3 , k p v = —ik p VT, (7) 

withX 3 ee {7 2 2 + 27 32 (27+r 13 +r 43 )|fi c | 2 /[7(7+r 3 )]} 1 / 2 . By taking a contour consisting 
of real axis and a semi-circle in the lower half complex plane [see the curves with arrows 
shown in figure [2]( a)], we can calculate the integration in equation ([6]) analytically by 
just calculating the residues corresponding to the two poles, and obtain exact result 
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Figure 2. (Color online) (a): Two poles (A2,— 1X3), (0, — ik v vx) of the integrand 
in equation ([6]) in the lower half complex plane. The closed curve with arrows is the 
contour chosen for calculating the integration in equation (j6|) by using residue theorem, 
(b): Absorption spectrum lm(K) as a function of w for the hot molecular system. The 
solid (dashed) line for \Cl c \ = 800 MHz (\Cl c \ = 0). Definitions of Im(K) min , Im(A") max , 
and the width of transparency window Ttw are indicated in the figure. 



for the integration. Since the expression is lengthy, we just write down the one with 
A 2 = A 3 = 0, Alo d > j jh 7: 

K = ^ + /Ci + /C 2 , (8) 

^n ls Au D {2-f 32 (u + i7 2 i)A(-iJ*f 3 ) - iX 3 B] 

1 " 7 r 37 (Ac4 - X|)X 3 [|n c | 2 - (oj + i 721 )(o; + 2X3)] ' 

y^i3 [2 732 (a; + z72i)A(-2Aa; D ) - iAc^p-B] 

2 ~ 7r 37 (^! - Au,f,)[|ft c | 2 - (w + i 721 )(o; + *Au; D )] ' 

where Acuq = fc p i>T (Doppler width), y4(/c p w) = X^'p— 7 |f2 c | 2 cx 22 and £> = 7 r 37 |f2 c | 2 <7 22 . 
Note that /Ci (/C 2 ) is contributed by the first (second) pole. For cold molecules the second 
pole in equation (J7|) does not exist, thus /C 2 = 0. However, for hot molecules one has 
/C 2 7^ due to Doppler effect, and hence the system may have very different quantum 
interference characters comparing with that of cold molecules. 

In most cases, K(w) can be Taylor expanded around the center frequency of the 
probe field (corresponding to w = 0), i.e., K(uS) = K + K\io + (l/2)K 2 u 2 . . ., where 
Kj = (&>K/ ' 800^)^=0. The coefficients K describes the phase shift (real part) and the 
absorption (imaginary part) per unit length and l/Ke(Ki) and l/Ke(K 2 ) represent the 
group velocity v g and group-velocity dispersion, respectively. 

3.2.1. Transparency window of probe-field absorption spectrum. Shown in figure [2(b) 
is Im(i^) as a function of 00. The dashed (solid) line is for \tt c \ = (|fi c | = 800 MHz). 
System parameters given by T 13 = T 23 = r 43 = 1.77 x 10 7 s" 1 , 7 = 0.47 x 10 6 s _1 , 
j]f = 4 x 10 6 s" 1 , Alu d = 1-22 GHz, « 13 = 5 x 10 10 cm^s" 1 and a£ = ovg = 0.5. 
One sees that the absorption spectrum of the probe field for \Q C \ = has only a single 
absorption peak. However, a transparency window opens for a |fi c = 800 MHz. The 
minimum (Im(i^) min ), maximum (lm(K) max ), and width of transparency window (r TW ) 
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are defined in the figure. 

From equation (jSJ), we obtain the minimum of lm(K) at u = 0: 

Aw D Vl + Xi 1 + Xtl + y/xJ 



where x = |f2 c | 2 7 31 / (jAcj^) and x x = |fi c | 2 / (72iAcjd) ar e two dimensionless parameters. 
It is interesting that the system has absorption and gain, reflected by the first and the 
second terms on the right hand side of equation ([9]) . The gain is due to non- vanishing 
7 and cr^L Obviously, if x 3> 1 and x\ ^> 1, i.e. |fi c | 2 73i 3> jAcj^j and |f2 c | 2 ^> 721 Au;.d, 
one has Im(J^) min m 0, i.e. a large and deep transparency widow in the absorption 
spectrum is opened. The inequalities can be taken as the EIT condition [T6J [27] of the 
system. When 721 ~ 7, this condition is simplified to |fi c | 2 73i ^ 7Ao;fj. 

Under the above condition, we obtain Im(i^) max ~ K^ally/n/ Au^ located at 
uj m ±f2 c , and 

1/2 

(10) 



3.2.2. EIT- ATS crossover. One of our main purposes is to explicitly analyze the 
detailed characters of quantum interference effect of the system, which can be done 
by extending the spectrum decomposition method introduced in |U EJ EJ [7J- Note that 
JCj (j = 1, 2) in equation (jHJ) can be decomposed as 

K^vJ^t + ^j-), (11) 
\uj-5 j+ u)-5j-J 

where rjj, Aj± are constants, 5j + and <5j-_ are two spectrum poles, all of which have 



been given explicitly in Appendix A From equation ( II ip we can get explicit expressions 
of Im(JCj) (j = 1,2). However, their general expressions are lengthy and complicated. 
In order to illustrate the quantum interference effect in a simple and clear way, we 
decompose Im(/C,-) according to different regions of Q c . 

(i). Weak control field region (i.e. |Q C | < f2 rc f = Aw/j/2): In this region, one has 
Re(<5j±)=0, Im(A,--i-)=0, we obtain 



2 2 

Im(A-) = 5>(C 3 ) = J> (-j^ + 

.7=1 .7 = 1 ■ 7 + 



L0 2 + Wf 



L t + L 2 , (12) 



where L\ and L2 are defined by 



T V2C: 



u 2 + W?_ u 2 + Wl 
V1C1+ V2C2+ 
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Figure 3. (Color online) EIT-ATS crossover for hot molecules in the co-propagating 
configuration, (a): Absorption spectrum in the region |fi c | < £7 rc f = 
contributed by positive L\ (dashed-dotted line), negative Li (dashed line), and total 
absorption spectrum Im(K) (solid line), (b) and (c): Absorption spectrum by two 
Lorentzians (dashed-dotted line), destructive interference (dashed line), and total 
absorption spectrum lm(K) (solid line), in the region |O c | > f2 ro f and |f2 c | f2 re f, 
respectively. Panels (a), (b) and (c) correspond to EIT, EIT-ATS crossover, and ATS, 
respectively. 



with real constants 



C j+ -- 
Cj-- 

w x ± 
w 2 ± 

1 1 — 



.(W j+ + TJ)/(W j 



i 



TJ)/(W j+ - W j 



1 2 



2 

721 - 
721 - 



X 3 + 72i ± V(X 3 -72i) 2 -4|fi c | 

Alo d + 72i ± y/{Au D -721) 2 -4|fi c | 2 
X 3 B 



2 l32 A(-iX 3 ) 
Au D B 



(14) 

2 l32 A(-iAcu D ) v ; 

Shown in figure [3]^a) are results of L±, which is a positive single peak (the dashed- 
dotted line), and L 2 , which is a negative single peak (the dashed line). System 

1.77 x 10 7 s-\ 7 = 0.47 x 10 6 s -1 , 
= 414 MHz. The sum of the positive 



parameters are given by T 



13 



23 



43 



= 4 x 10 e 



Au 



D 



1.22 GHz, and fi„ 



L\ and negative L 2 gives lm(K) (the solid line), which displays a absorption doublet 
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with a significant transparency window near at u — 0. Because there exists a destructive 
interference in the probe-field absorption spectrum, the phenomenon found here belongs 
to EIT according to the criterion given in [HJ EJ [7] . 

(ii). Intermediate control field region (i.e. |f2 c | > fi re f) : By extending the approach 
by Agarwal [I], we can decompose Im(/Cj) (j — 1, 2) as 



Im(/Cj 



[CO 



<5p 2 + W] 



+ 



25^ L 



5p 2 



w + 



iy 2 



(15) 



where 



Wi 

w 2 

61 = 
5 r 2 = 

91 = 

92 = 



( 721 + X 3 )/2, 
(72i + Au D )/2, 



4|fi c | 2 -(X 3 - 7 2i) 2 /2, 



-^3 - 721 



(Au D - 7 2i) 2 /2, 
X 3 B 



2 2 732 A(-^X 3 ) 
Aw D - 72i , Au D B 



2 ' 2-y sa A(-iAu) D y (16) 

The first two terms in the first square bracket on the right hand side of equation (115jl are 
two Lorentzians, resulted from the absorption from two different pathways corresponding 
to the two dressed states created by the coupling field. The terms in the second square 
bracket are interference terms, the magnitudes of which are controlled by the parameter 
gj. If gj > (<7j < 0) the interference is destructive (constructive). 

Figure E^b) shows the result of the probe-field absorption spectrum as functions 
of u for \Q C \ > Q ve f. The dashed-dotted line (dashed line) denotes the contribution by 
two Lorentzians (interference terms). We see that the interference is destructive. The 
solid line gives the result of Im(i^). System parameters used are the same as those in 
panel (a) but with il c — 1 GHz. A transparency window opens due to the combined 
effect of EIT and ATS, which is deeper and wider than that in panel (a). We call such 
phenomenon as EIT-ATS crossover. 

(iii). Large control field region (i.e. \Q C \ ^> f2 rc f): In this case, the quantum 
interference strength gj/Sj in equation (fl5|) is very weak and negligible. We have 



Im(/C,) « | 



Wn 



Wj 



[CO 



5p 2 + W? 



(17) 



(cu + 5p 2 + W}\ ' 

being to a sum of two Lorentzians. 

Shown in the panel (c) of figure[3]is the result of the probe-field absorption spectrum 
as functions of u for |Q C | ^> f2 rc f. The dashed-dotted line represents the contribution by 
the sum of the two Lorentzians. For illustration, we have also plotted the contribution 
from the small interference terms [neglected in equation (flTI) ]. denoted by the dashed 
line. We see that the interference is still destructive but very small. The solid line is the 
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1 2 3 4 5 



1 c 1 ref 

Figure 4. (Color online) The "phase diagram" illustrating the transition from 
EIT to ATS for hot molecules in the co-propagating configuration. Shown is 
Jm(K) u=0 /lm(K) max as a function of |fi c |/fi ref . Three regions (EIT, EIT-ATS 
crossover, and ATS) are divided by two dash-dotted lines. 

curve of Im(A'), which has two resonances at u « ±Q C . Parameters used are the same 
as those in panel (a) and (b) but with Q c = 3 GHz. Obviously, the phenomenon found 
in this situation belongs to ATS because the transparency window opened is mainly due 
to the contribution of the two Lorenztians. 

From above results, we see that the probe-field absorption spectrum experiences 
a transition from EIT to ATS as Q c is changed from weak to strong values. Since in 
three-level systems such phenomenon happens quite often and is universal, we divide 
quantum interference effects into three classes, i.e. the EIT region (|fi c | < £l T ef), the 
region of the EIT-ATS crossover (1 < |f2 c |/f2 rcf < 4), and ATS region (|O c |/O re f > 4). 
Figure |4] shows a "phase diagram" that illustrates the transition from the EIT to ATS 
by plotting Im(A) w= o/Im(A) max as a function of |fi c |/f2 re f. Note that we have defined 
Im(A) aJ=0 /Im(A) max = 0.01 as the border between EIT-ATS crossover and ATS regions. 

3.2.3. Comparison with experiment. To check the theoretical prediction given above, 
it is necessary to make a comparison with the experiment reported recently by Lazoudis 
et al. [9], which was performed with a co-propagating configuration. Using system 
parameters = r 23 = r 43 = 1.77 x 10 7 s -1 , 7 = 0.47 MHz, 7^ = 4 MHz, and 
Acu£> = 1.22 GHz, we have calculated probe-field absorption spectrum Im(A) as a 
function of frequency u, with Q c = 414 MHz (EIT region) and the control-field detuning 
—55 MHz. The result is plotted as the dashed line of figure El which agrees fairly with the 
experimental one (the solid line) measured in [9] (see figure 5(a) of [9]). Note that here 
we have plotted the quantity Im(A'), which is proportional to fluorescence intensity 
(measured in [9]) related to the state |3) because 033 ~ 2|f2 p | 2 Im(A)/(7 + T 3 ) [29J. 
The small difference for depth and width of the EIT dip between our result and the 
experiment is due to the approximation by using the modified Lorentzian distribution 
to replace the Maxwellian velocity distribution. 
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Figure 5. (Color online) Probe-field absorption spectrum Im(i ; i')/Im(A') max as a 
function of frequency u>, with f2 c = 414 MHz (EIT region). The dashed line is 
theoretical result. The solid line is the experimental one reported in Ref. [pj. 



3.3. Hot molecules with counter-propagating configuration 

We now move to the situation when the probe and control fields are arranged as a 
counter-propagating configuration. Here, d 2 \ = A 2 — A x — 2k p v + 2721 and d 32 = 
A3 — A2 + k p v + «732- Then we obtain 

c 



K = - 



7r 37 



(£1 + /c 2 



£1 



/c 2 



^^[2732(0; + z2X 3 )A(-zX 3 ) + zX 3 £?] 
(Aoo 2 D - Xf)X 3 [|fi c | 2 - (w + 2iX 3 )(ou + iX 3 )] ' 

^[2732(0; + i2Auj D )A(-iAcu D ) + zAu^-B] 
(XI - AcuiypJ 2 - (w + i2Auj d )(uj + iAu D )} 

A 3 - iX 3 



and 



where /Ci and /C2 are obtained from the poles k p v 
respectively. 

We have carried out a similar spectrum decomposition as that did for the co- 
propagating configuration given above. For saving space, here we omit concrete 
expressions of the spectrum decomposition but present probe-field absorption spectra 
in three typical control-field regions in figure [6j 

Shown in the panel (a) of figure [6] is the result of probe-field absorption spectrum 
Im(i^) in weak control-field region (i.e. \Q C \ < fi re f) as a function of a; for Q c = 500 MHz. 
As in figure EJ^a), Im(J^) is also the sum of two terms, i.e. L\ and L 2 . Nevertheless, now 
both L\ and L 2 are positive, as illustrated by the dashed-dotted line and dashed line, 
respectively. We see that lm(K) (the solid line) displays only a positive single peak, there 
is no transparency window, and the reason is that the quantum interference becomes 
constructive (the red dashed line) for the counter-propagating configuration. Thus, 
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Figure 6. (Color online) Probe-field absorption spectrum for hot molecules in the 
counter-propagating configuration, (a): Absorption spectrum in the region |f2 c | < 
57 rc f = Acjd/2 contributed by positive L\ (dashed-dotted line) and L2 (dashed line), 
and total absorption spectrum Im(A") (solid line), (b) and (c): Absorption spectrum 
by two Lorentzians (dashed-dotted line), constructive interference (dashed line), and 
total absorption spectrum Im(A') (solid line), in the region |fi c | > fi re f and |£l c | 3> f2 ra f , 
respectively. 



different from the case of the co-propagating configuration, in weak control-field region 
an EIT which we have defined as transparency window plus a destructive interfrence 
does not exists. 

Shown in figure [6](b) and (c) are results of the probe-field absorption spectra as 
functions of u for \Q C \ > f2 rc f and |O c | ^> f2 rc f, respectively. System parameters are 
given by T 13 = T 23 = T 43 = 1.77 x 10 7 s _1 , 7 = 0.47 x 10 6 s~\ 7™ 1 = 4 x 10 6 s -1 , 
and Aud = 1.22 GHz, with Q c = 3 GHz (in the intermediate control-field region) and 
Q c = 25 GHz (in the large control-field region) for the panel (b) and the panel (c), 
respectively. The dashed-dotted line (dashed line) denotes the contribution by the sum 
of two Lorentzians terms (interference terms) in Im(J^). The solid line gives the result of 
Im(i\~). We see that the interferences near the probe- field center frequency (i.e. u = 0) 
are always constructive. Consequently, different from the case of the co-propagating 
configuration, no EIT- ATS crossover happens. 

Shown in figure [7] is the "phase diagram" that illustrates the transition from the 
constructive interference to ATS for the counter-propagating configuration by plotting 
Im(i^) w=0 /Im(i ; r) max as a function of |f2 c |/f2 ref . Three regions are divided as constructive 
interference, ATS with constructive interference, and ATS, respectively. 
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Figure 7. (Color online) lm(K) ul= o/Im(K)j nax as the function of the control field 
|r2 c |/f2 ro f for hot molecules in the counter-propagating configuration. Three regions 
(constructive interference, ATS with constructive interference and ATS) are divided 
by two dashed-dotted lines. 



3.4- Cold molecules and comparison for various cases 

Our model presented in section [2]is also valid for cold molecules. In this case, one should 
take v — in the Bloch equation (j2J), and f(v) = 5(v) in the Maxwell equation fl3]). 
The solutions (jl]) and fl5]) are still valid but one must take v — there. However, the 
dispersion relation (jfJJ) is replaced by 



with T = 721 + |f2 c | 2 ((733 — C22 )/[732(°"n — °33 )]■ Here A2 = A3 = has been taken 
for simplicity. 

A similar spectrum decomposition can be done like that did for hot molecules, 
which is omitted here. Shown in figure M is the probe-field absorption spectrum 
Im(i\~) w=u /Im(i\~) max as a function of the control field |f2 c |/f2 re f, where f2 rc f = |^y 2 i — 
731 1/2. System parameters are given by Ti 3 = T 2 3 = r 43 = 1.77 x 10 7 s _1 , 7^°' = 1 x 10 3 
s -1 and all = F rom the figure, we obtain the similar conclusion as that obtained for 
co-propagation configuration, i.e. the probe-field absorption spectrum experiences also 
a transition from EIT to ATS as Q c is increased. The quantum interference effect in 
the system can divided into three regions, i.e. the EIT region (|O c | < f2 re f), the region 
of the EIT-ATS crossover (1 < |f2 c |/f2 re f < 5), and ATS region (|f2 c |/Q re f > 5). 

From the results given above, we see that the quantum coherence in the open A- 
type molecular system has very interesting features, depending on the existence or non- 
existence of the Doppler broadening, and also depending on the beam propagating (co- 
propagating or counter-propagating) configurations. For comparison, in Table [1] some 
useful physical quantities, including EIT condition, absorption spectrum lm(K)\ w= o, 
group velocity v g , and width of transparency window Ttw? are presented for several 
different physical cases. 




(19) 
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Figure 8. (Color online) Transition from EIT to ATS for cold molecules. Shown is 
Im(A') w= o/Im(A" ) max as a function of |f2 c |/f2 re f, where il rc f = \j2i — 73i|/2- Three 
regions (EIT, EIT- ATS crossover and ATS) are divided by two dash-dotted lines. 



Table 1. Propagating properties of the probe held for various open A- type 
molecular systems, including EIT condition, absorption spectrum Im(.K^)| w= o, width 
of transparency window Ttw, and group velocity v g for three different cases. Other 
quantities appeared in the Table have been defined in the text. Mol.=Molcculcs, Co- 
prop. =Co-propagating configuration, Cou.-prop.=Counter-propagating configuration. 



System EIT condition Im(iT)| w=0 T TW v 



9 



Hot Mol. (Co-prop.) ^ < \Q C \ 2 < ^™ « 

Hot Mol. (Cou.-prop.) no EIT ^wP 2 2\Q C \ - Au D 



Cold Mol. 72i73i < \tt c \ 2 < 



\n c \ 2 1 c| ^ u v^«i3 

7|i ^13721 2|r2 e | 2 |a£ 

4 \fl c \ 2 731 «M3 



The first line in the Table is for hot molecules working in the co-propagating 
configuration; the second line is for hot molecules working in the counter-propagating 
configuration; the third line is for cold molecules. There are EIT, EIT-ATS crossover, 
and ATS for both cold molecules and the hot molecules with the co-propagating 
configuration. But there is no EIT and no EIT-ATS crossover for the hot molecules 
with the counter-propagating configuration. Experimentally, up to now only the EIT in 
the co-propagating configuration has been demonstrated recently by experiment [9]. 

4. Nonlinear pulse propagation 

The theoretical approach given in the last two sections is valid not only for continuous- 
wave but also for pulsed probe fields. However, if the probe field is pulsed and has a 
larger amplitude, nonlinear effect induced by Kerr nonlinearity inherent in the system 
must taken into account. We stress that the theoretical scheme proposed in the present 
work is very suitable for the study of pulse propagation in multi-level systems. 

In this section, we investigate nonlinear pulse propagation, especially ultraslow 
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optical solitons, in the present open hot molecular system with co-propagating 
configuration by using the method of multiple-scales. For this aim, we take the 
asymptotic expansion crfl-crfi = Y^ m =i,2,- eTna jT\ = E m =i,2 r . e m Q p m \ with crvV = 
and 032^ = 0, where e is a small parameter denoting the typical amplitude of f2 p and 
all quantities on the right hand side of the asymptotic expansion are considered as 
functions of the multi-scale variables z m = e m z (m = 0,1,2), t m = e m t (m = 0,1). 
Substituting the expansion into the MB equations ([2]) and ([3]), we obtain a series of 
linear but inhomogeneous equations for and Vt p m ^ (m = 1-4), which can be solved 
order by order. 

The zeroth-order (m = 0) and the first-order (m = 1) solutions are the same 
as that given respectively by equation (J4]) and ([5]), by now 9 = K(u)z — ut and F 
is yet to be determined envelope function of the "slow" variables t\, z\ and z 2 - In 
the second order (m = 2), a divergence- free solution for flp requires the solvability 
condition i[dF/dzi + (dK/du)dF/dti\ = 0, which shows that the envelope function 
F travels with complex group velocity (dK/du)" 1 . Explicit expressions of the second 



order solution have been given in Appendix B 



In the third order (m = 3), the Kerr nonlinearity of the system plays a role. A 

(3) 

divergence-free solution for Q p gives rise to the equation 
dF ld 2 Kd 2 F 



i- 



'dz 2 2 dco 2 dt\ 
where a = lm(K) = e 2 a and 



W\F\ 2 Fe~ mz2 = 0, (20) 



w = _ ku r> m ^ + f+^f-f , ( 21) 

J-oo \U C \ 2 - (u + d 2 i){uj + d 3 i) 

,(2) .(2) nTW , J2) 



with coefficients a u , a 22 and defined in|Appendix B 



Combining equation (l2"0|) and the solvability condition in the second order, we 
obtain 

<9 1 8 2 K B 2 JJ 

where r = t—z/v g and U = eF. Equation (|22p is a nonlinear Schrodinger (NLS) equation 
describing time evolution of the envelope function F, in which W is proportional to 
third-order nonlinear susceptibility (Kerr coefficient) relevant to self-phase modulation, 
which is necessary for the formation of a shape-preserved probe pulse. 

The key for the formation and propagation of an optical soliton in the system 
requires two conditions. The first is a balance between dispersion and nonlinearity, and 
the second is the absorption of the probe field must be negligibly small. Generally, 
the coefficients of the equation (1221) are complex, which means that a soliton, even 
if it is produced initially, may be highly unstable during propagation. However, as 
shown below, a realistic set of system parameters can be found under the EIT condition 
so that the imaginary part of these coefficients can be made much smaller than their 
corresponding real part. Thus it is possible to get a shape-preserving nonlinear localized 
solution that can propagate a rather long distance without a significant distortion. 
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Figure 9. (Color online) (a): Ultraslow optical solitons and their interaction in the 
hot molecular system, (a): Three-dimensional plot of the waveshape of |O p /[/o| 2 as a 
function of z/Ld and i/ r o- (b) : Collision between two ultraslow optical solitons. 



Neglecting the small imaginary part of the coefficients and taking to — 0, equation 
( 1221 can be written into the dimensionless form idu/ds + d 2 u/da 2 + 2\u\ 2 u = 0, with 
s = —z/(2L D ), a = t/t , and u = U/U . Here r is typical pulse duration, L D = Tq/K 2 



is typical dispersion length, and Uq = (1/V ) y K 2 /W is typical half Rabi frequency of 

the probe field, with K 2 and W being the real part of K 2 = (d 2 K/du 2 ) w= o and ff^^o, 
respectively. Then one can obtain various soliton solutions for u. A single-soliton 
solution in terms of the half Rabi frequency reads 



with Kq = Re(if)| w= o, which describes a bright soliton traveling with the propagating 
velocity v g = [Re(dK /du)}' 1 ]^^. 

We now give a realistic parameter set for the formation of the optical soliton given 
above. For a hot Li 2 molecular gas, we choose Q c = 600 MHz, A 2 = A 3 ps 2.36 x 10 7 
s" 1 , r = 1.0 x 10 -7 s, u p = 4.46 x 10 14 s" 1 , and other parameters are the same as those 
given in the previous text. Then we obtain K 2 = (5.51 + 0.672z) x 10~ 16 cm~ 1 s 2 and 
W = (1.75 + 0.298z) x 10~ 16 chi'V, L d = L NL = 18.2 cm, and U = 1.77 x 10 7 s" 1 . 
One sees that the imaginary part of K 2 and W is indeed much smaller than their 
corresponding real part. The reason of so small imaginary part is due to the quantum 
interference effect contributed by the control field. 

The propagating velocity of the probe pulse can be estimated by the real part of the 
linear dispersion relation ([6]). At the probe- field center frequency (i.e. u = 0) we obtain 
v g = \Re(dK/duj)\ ll j = o\~ 1 ~ 2.13 x 10~ 4 c. Consequently, the optical soliton obtained 
may travel with an ultraslow propagating velocity in the system. 

The stability of the ultraslow optical soliton described above can be checked by using 
numerical simulations. In figure E{a), we show the wave shape of \Q p /Uq\ 2 as a function 
of z/Ld and t/ro. The solution is obtained by numerically solving Eq. ( 1221) with full 
complex coefficients included. The initial condition is given by fl p (0,t) = U sech(t / r ) . 
We see that the amplitude of the soliton undergoes only a slight decrease and its width 





(23) 
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undergoes a slight increase due to the influence of the imaginary part of the coefficients. 
A simulation of the interaction between two ultraslow optical solitons is also carried 
out by inputting two identical solitons [see figure |9] (b)]. The initial condition is 
Q p (0,t) = t/ sech(t/r — 5) + U sech(t/T + 5). As time goes on, they collide, pass 
through, and depart from each other. The two solitons recover their initial waveforms 
after the collision. However, a phase shift is observed after the collision. 



We have developed a systematic analytical approach on linear and nonlinear pulse 
propagations in an open A-type molecular system with Doppler broadening. In linear 
case, by using residue theorem and spectrum decomposition method, we have proved 
that there exists a crossover from EIT to ATS for the co-propagating configuration. 
However, there is no EIT and hence no EIT-ATS crossover for the counter-propagating 
configuration. We have provided various explicit formulas, including probe-field 
spectrum decomposition, EIT condition, and width of EIT transparency window, as 
well as a comparison with the result of cold molecules. Our analytical result agrees well 
with the experimental one reported recently by Lazoudis et al [S]. In nonlinear case, 
by using the method of multiple-scales, we have derived a nonlinear envelope equation 
for probe-field propagation. We show that stable ultraslow solitons can be realized 
in the open molecular system. New theoretical predictions presented in this work 
are helpful for guiding new experimental findings in coherent molecular systems and 
may have promising practical applications in coherent molecular spectroscopy, precision 
measurement, molecular quantum state control, nonlinear pulse propagation, and so on. 
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Appendix A. Expressions of rjj, Aj±, and 5j± 



5. CONCLUSION 
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Appendix B. Second-order solution of MB Equations 
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(B.5) 

(B.6) 



with D = \Q C \ 2 - (u + d 2l )(u + d 31 ) and D x = 7 [( 7 + T 23 + r 43 )( 7 + T 31 ) + 
d 32 )(u + d* 32 ) + 2 732 [(2 7 + r 43 )( 7 + T 31 ) + 7 ri 3 ]|fi c | 2 . a£ } and 4? have been 
equation ([5]). 
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